Coulomb Parameter U and Correlation Strength in LaFeAsO 
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First principles constrained density functional theory scheme in Wannier functions formalism has 
been used to calculate Coulomb repulsion U and Hund's exchange J parameters for iron 3d electrons 
in LaFeAsO. Results strongly depend on the basis set used in calculations: when 0-2p, As-4p, and 
Fe-3d orbitals and corresponding bands are included, computation results in U =3-^4 eV, however, 
with the basis set restricted to Fe-3d orbitals and bands only, computation gives parameters corre- 
sponding to F°=0.8 eV, J=0.5 eV. LDA+DMFT (the Local Density Approximation combined with 
the Dynamical Mean-Field Theory) calculation with this parameters results in weakly correlated 
electronic structure that is in agreement with X-ray experimental spectra. 
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Following the discovery of high- Tc superconductivity in 
iron oxypnictide La0i_2;Fa;FeAs [l|, a question of the in- 
fluence of electronic correlation effects on the normal and 
superconducting properties of LaFeAsO has arisen. In 
striking similarity with high-T^ cuprates, undoped mate- 
rial LaFeAsO is not superconducting with antifcrromag- 
netic commensurate spin density wave developing below 
150 K y. Only when electrons (or holes) are added to the 
system via doping, antiferromagnctism is suppressed and 
superconductivity appears. As it is generally accepted 
that Coulomb correlations between copper 3d electrons 
are responsible for cuprates anomalous properties, it is 
tempting to suggest that the same is true for iron 3d 
electrons in LaFeAsO. 

Correlation strength in a system is determined by ra- 
tio of Coulomb interaction parameter U and band width 
W . If U /W is significantly less than 1 then the system is 
weakly correlated and results of the Density Functional 
Theory (DFT) calculations are reliable enough to explain 
its electronic and magnetic properties. However, if U 
value is comparable with W or even larger then the sys- 
tem is in intermediate or strongly correlated regime and 
Coulomb interactions must be explicitly treated in elec- 
tronic structure calculations. For LaFeAsO the bands 
formed by Fe-3(i states have width w4 eV (see shaded 
area in the lower panel of Fig. [1]) , so an estimation for 
Coulomb interaction parameter U should be compared 
with this value. 

In practical calculations, U is often considered as a 
free parameter to achieve the best agreement of calcu- 
lated and measured properties of investigated system. 
Sometimes U value could be estimated from the experi- 
mental spectra. The first principles justified methods to 
determine Coulomb interaction parameter U value are 
constrained DFT scheme [3|, where in DFT calculations 
the d-orbital occupancies are fixed to the certain values 
and U is numerically determined as a derivative of d- 



orbital energy over its occupancy and Random Phase 
Approximation (RPA) method, where screened Coulomb 
interaction between d-electrons is calculated via pertur- 
bation theory [4[. In Ref. |5| it was proposed to use 
in LaFeAsO U=A cV obtained in RPA calculations for 
metallic iron [6| . 

This value for Coulomb parameter (with Hund's ex- 
change parameter J=0.7 eV) was used in Dynami- 
cal Mean-Field Theory (DMFT) 0] calculations for 
LaFeAsO [5|, |8|, |9|. Results of these works show iron 
3d electrons being in intermediate or strongly correlated 
regime, as it is natural to be expected for Coulomb pa- 
rameter value t/=4 eV and Fe-3d band width «4 eV. 

The most direct way to estimate correlation effects 
strength in a system under consideration is to compare 
the experimental spectra with densities of states (DOS) 
obtained in DFT calculations. For strongly correlated 
materials additional features in the experimental pho- 
toemission and absorption spectra appear that are inter- 
preted as lower and upper Hubbard bands absent in the 
DFT DOS. If no such additional features are observed 
and DOS obtained in DFT calculations satisfactorily de- 
scribe the experimental spectra then the material is in 
weakly correlated regime. 

LaFeAsO was studied by soft X-ray absorption and 
emission spectroscopy 10], X-ray absorption spec- 
troscopy (O K-edge) _11|] and photoemission spec- 
troscopy I12|. In all these works the conclusion was that 
DOS obtained in DFT calculations gave good agreement 
with the experimental spectra and the estimations for 
Coulomb parameter value are U <\ eV ll|. That con- 
tradiction with results of the DMFT calculations [a, la ISl 
shows that first principles calculation of Coulomb inter- 
action parameter U value for LaFeAsO is needed to de- 
termine the correlation effects strength in this material. 
Results of such calculations by constrained DFT calcula- 
tions are reported in the present work. We have obtained 



the value U <1 eV for Fe-Sd band that agrees with the 
estimates from spectroscopy. Recently the RPA calcula- 
tions for Coulomb interaction parameter U in LaFeAsO 
were reported, where U was estimated as 1.8-^2.7 eV [l3| . 
It is important to note that Coulomb interaction pa- 
rameter U value depends on the model where it will be 
used and, more precisely, on the choice of the orbital 
set that is taken explicitly into account in the model. 
For example, in constrained DFT calculations for high- Tc 
cuprates the resulting U value for Cu d-shell was found 
between 8 and 10 eV 1J|. The U value in this range 
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was used in cluster calculations where all Cu d-orbitals 
and p-orbitals of neighboring oxygens were taken into ac- 
count and calculated spectra agree well with experimen- 
tal data [15| . However, in one band model, where only 
x^ — y^ orbital per cooper atom is explicitly included in 
the calculations, the U value giving good agreement with 
experimental data falls down to 2.5-^3.6 eV [l6|, that is 
3-4 times smaller than constrained DFT value. 

The same situation occurs for titanium and vanadium 
oxides: the U value from constrained DFT calculations 
is «6 eV and cluster calculations where all d-orbitals and 
p-orbitals of neighboring oxygens were taken into account 
with U close to this value gave good agreement between 
calculated and experimental spectra [l7[. However, in 
the model where only partially filled t2g orbitals are in- 
cluded, much smaller U value (corresponding to Slater 
integral F'^—S.b eV) gives the results in agreement with 
experimental data [l8i | . 

It is interesting that such a small U value can be ob- 
tained in constrained DFT calculations for titanates and 
vanadates where only i2g-orbital occupancies are fixed 
while all the other states (cg-orbitals of vanadium and 
p-orbitals of oxygens) allowed to relax in self-consistent 
iterations [isl, ll9|. So the calculation scheme used in 
constrained DFT (the set of the orbitals with fixed occu- 
pancies) should be consistent with basis set of the model 
where the calculated U value will be used. 

Another source of uncertainty in constrained DFT 
calculation scheme is a definition of atomic orbitals 
whose occupancies are fixed and energy calculated. In 
some DFT methods, like Linearized MufRn-Tin Orbitals 
(LMTO) ^, these orbitals could be identified with 
LMTO. However, in other DFT calculation schemes, 
where plane waves are used as a basis, like in pseudopo- 
tential method [2l'| one should use more general definition 
for localized atomic like orbitals such as Wannier func- 
tions (WFs) [231 . The practical way to calculate WFs for 
specific materials using projection of atomic orbitals on 
Bloch functions was developed in [23|. 

In Fig. [1] the total and partial DOS for LaFeAsO ob- 
tained in LMTO calculations are shown. Crystal field 
splitting for Fe-3d orbitals in this material is rather weak 
(Ac/=0.25 eV) and all five d orbitals of iron form com- 
mon band in the energy region (—2, -1-2) eV relative to 
the Fermi level (see grey region on the bottom panel in 
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FIG. 1: Total and partial densities of states for LaFeAsO 
obtained in DFT calculation in frame of LMTO method. 



Fig. [T]) . There is a strong hybridization of iron t2g orbitals 
with p orbitals of arsenic atoms which form nearest neigh- 
bors tetrahedron around iron ion. This effect becomes 
apparent in the energy interval (—3, —2) eV (white re- 
gion on the bottom panel in Fig. [T]) where band formed 
by p orbitals of arsenic is situated. More week hybridiza- 
tion with oxygen p states reveals in (—5.5, —3) eV energy 
window (black region on the bottom panel in Fig. [T]). 

We have calculated Coulomb interaction U and Hund's 
exchange J parameters for WFs basis set via constrained 
DFT procedure with fixed occupancies for WFs of d sym- 
metry. For this purpose we have used two calculation 
schemes based on the different DFT methods. One of 
them involves linearized muffin-tin orbitals produced by 

the TB-LMTO-ASA code H; corresponding WFs cal- 

I i 
culation procedure is described in details in ReL 24. The 

second one is based on the plane waves obtained within 
the pseudopotential plane-wave method PWSCF, as im- 
plemented in the Quantum ESPRESSO package [21[, and 
described in details in Ref. |25|. The difference between 
the results of these two schemes could give an estimation 
for the error of U and J determination. 

The WFs are defined by the choice of Bloch functions 
Hilbert space and by a set of trial localized orbitals that 
will be projected on these Bloch functions. We performed 
calculations for two different choices of Bloch functions 
and atomic orbitals. One of them includes only bands 
predominantly formed by Fe-3d orbitals in the energy 
window (—2, +2) eV and equal number of Fe-Sd orbitals 
to be projected on the Bloch functions for these bands. 
That choice corresponds to the model where only five d- 
orbital per Fe site are included but all arsenic and oxygen 
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FIG. 2: (Color online) Module square of d^2_y2-like Wannier 
function computed for Fe-3d bands only (left panel) and for 
full set of 0-2p, As-4p and Fe-Sd bands (right panel). Big 
sphere in the center marks Fe ion position and four small 
spheres around it correspond to As neighbors. 



p-orbitals are omitted. Second choice includes all bands 
in energy window (—5.5, +2) eV that are formed by O- 
2p, As-4p and Fe-Sd states and correspondingly full set of 
0-2p, As-4p and Fe-3d atomic orbitals to be projected on 
Bloch functions for these bands. That would correspond 
to the extended model where in addition to d-orbitals all 
p-orbitals are included too. 

In both cases we obtained Hamiltonian in WF basis 
that reproduces exactly bands predominantly formed by 
Fe~3d states in the energy window (—2, +2) eV (Fig. [T]), 
but in the second case in addition to that bands formed 
by p-orbitals in the energy window (—5.5, —2) eV will be 
reproduced too. However, WFs with d-orbital symmetry 
computed in those two cases have very different spatial 
distribution. In Fig. [2] the module square of dx^-y^-^i^e 
WF is plotted. While for the case when full set of bands 
and atomic orbitals was used (right panel) WF is nearly 
pure atomic d-orbital (iron states contribute 99%), WF 
computed using Fe-Sd bands only is much more extended 
in space (left panel). It has significant weight on neigh- 
boring As ions with only 67% contribution from central 
iron atom. 

The physical reason for such effect is p — d hybridiza- 
tion that is treated explicitly in the case where both p- 
and d-orbitals are included. In the case where only Fe-3d 
bands are included in calculation p — d hybridization re- 
veals itself in the shape of WF. Fe-3d bands in the energy 
window (—2, +2) eV correspond to antibonding combi- 
nation of Fe-3(i and As-4p states and that is clearly seen 
on the left panel of Fig. [21 

The different spatial distribution for two WFs calcu- 
lated with full and restricted orbital bases can be ex- 
pected to lead to different effective Coulomb interaction 
for electrons occupying these states. The results of con- 
strained DFT calculations of the average Coulomb inter- 
action U and Hund's exchange J parameters for electrons 
on WFs computed with two different set of bands and 
orbitals (and using two different DFT methods: LMTO 
and pseudopotential) are presented in Tab. [T] 



TABLE I: The constrained DFT calculated values of average 
Coulomb interaction U and Hund's exchange J (eV) param- 
eters for d-symmetry Wannier functions computed with two 
different sets of bands and orbitals. 



DFT method 



TB-LMTO-ASA 
PWSCF 



separate Fe-3d band 



(7=0.49, J=0.51 
(7=0.59, J=0.53 



full bands set 



U=3.10, J=0.81 
(7=4.00, J=1.02 



One can see that very different Coulomb interaction 
strength is obtained for separate Fe-Sd band and full 
bands set calculations. While the latter gives value 
3-^4 eV, that agrees with previously used values [a, |8|, |9|] , 
separate Fe-3(i band calculation results in 0.5-^0.6 eV, 
that is much smaller but agrees with spectroscopy esti- 
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mations 

The main reason for such a drastic difference between 
two calculations is very different spatial extension of the 
two WFs (see Fig. [2]): nearly complete localization on 
central iron atom for "full bands set" WF (99%) and only 
67% for "Fe-3d band set" WF. Another possible reason 
for strong reduction of the calculated U value in going 
from "full bands set" to "Fe-3(i band set" WF is screening 
via p — d hybridization with As-4p band that is situated 
just below Fe-3d band (see Fig.[T]). The effect of decreas- 
ing of the effective U value in several times going from 
full orbital model to restricted basis was found previously 
for high-Te cuprates (C/=8-^I0 eV for full p — d-orbitals 
basis 14 1 and 2.5-^3.6 eV for one-band model (l6J|). 

In constrained DFT calculations one obtains an av- 
erage Coulomb interaction U that can be estimated as 
U = F^ — J /2. Hence, Slater integral F" can be calcu- 
lated &sF° = U + J/2%. For "Fe-3d band set" WF that 
gives F°=0.8 eV at J=0.5 eV. With this set of parame- 
ters we performed the LDA-I-DMFT [20| calculations (for 
detailed description of the present computation scheme 
see Ref. l24h . The DFT band structure was calculated 
within the TB-LMTO-ASA method 0. Crystal struc- 
ture parameters were taken from Ref. 111. 

The restricted basis set including only Fe-3(i WFs 
was used in the LDA-fDMFT calculations. The effec- 
tive impurity model for the DMFT was solved by the 
QMC method in Hirsh-Fye algorithm [23]. Calcula- 
tions were performed for the value of inverse temperature 
(3—\Q eV~^. Inverse temperature interval < r < /3 was 
divided into 100 slices. 6 • lO*' QMC sweeps were used 
in self-consistency loop within the LDA+DMFT scheme 
and 12 ■ 10^ of QMC sweeps were used to calculate the 
spectral functions. 

The iron 3d orbitally resolved spectral functions ob- 
tained within DFT and LDA+DMFT calculations are 
presented in Fig. [3l The influence of correlation effects 
on the electronic structure of LaFeAsO is minimal: there 
are small changes of peak positions for 3z^ — t^ , xy and 
x^ — y^ orbitals (the shift toward the Fermi energy) and 
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FIG. 3: Partial densities of states for different Fe-3d orbitals 
obtained within the DFT (fiUed areas) and LDA+DMFT or- 
bitaUy resolved spectral functions for ^"=0.8 eV, J— 0.5 eV 
(bold lines). 



practically unchanged picture of spectral function dis- 
tribution for yz, zx bands. There are no appearance of 
neither quasiparticle peak on the Fermi level nor Hub- 
bard bands in the energy spectrum with such values of 
U and J. Hence LaFeAsO can be considered as weakly 
correlated material. 

This agrees with the results of soft X-ray absorption 
and emission spectroscopy study !lO!|. It was concluded 
there that LaFeAsO does not represent strongly corre- 
lated system since Fe L3 X-ray emission spectra do not 
show any features that would indicate the presence of the 
low Hubbard band or the quasiparticlepeak that were 
predicted by the LDA-FDMFT analysis [i,[i,i| with the 
large 11=4 eV. A comparison of the X-ray absorption 
spectra (O K-edge) with the LDA calculatioiis gave an 
upper limit of the on-site Hubbard [/ wl eV [11|. Pho- 



toemission spectroscopy study of LaFeAsO suggests [12 1 
that the line shapes of Fe 2p core-level spectra corre- 
spond to an itinerant character of Fe 3d electrons. It 
was demonstrated there that the valence-band spectra 
are generally consistent with band-structure calculations 
except for the shifts of Fe Sd-derived peaks toward the 
Fermi level. 

In conclusion, we have calculated the values of U and J 
via constrained DFT procedure in the basis of WFs. For 
minimal model including only Fe-3ci orbitals we have ob- 
tained Coulomb parameters F^=0.8 eV, J=0.5 eV. The 



LDA-fDMFT calculation with these parameters results 
in weakly correlated nature of iron d bands in this com- 
pound. This conclusion is supported by several X-ray 
spectroscopic investigations of this material. 
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